// Software License Agreement (BSD License)
//
// Copyright (c) 2012-2013, Fraunhofer FKIE/US
// All rights reserved.
// Author: Torsten Fiolka
//
// Redistribution and use in source and binary forms, with or without
// modification, are permitted provided that the following conditions
// are met:
//
//  * Redistributions of source code must retain the above copyright
//    notice, this list of conditions and the following disclaimer.
//  * Redistributions in binary form must reproduce the above
//    copyright notice, this list of conditions and the following
//    disclaimer in the documentation and/or other materials provided
//    with the distribution.
//  * Neither the name of Fraunhofer FKIE nor the names of its
//    contributors may be used to endorse or promote products derived
//    from this software without specific prior written permission.
//
// THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
// "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
// LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
// FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
// COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
// INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
// BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
// LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
// CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
// LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
// ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
// POSSIBILITY OF SUCH DAMAGE.

#include <include/sure/normal/base_histogram.h>

std::vector<sure::NormalType, sure::NormalTypeAllocator> sure::normal::initVectors(unsigned numberOfPolarBins)
{
  std::vector<NormalType, NormalTypeAllocator> vecs;
  const float polarBinSize = M_PI / float(numberOfPolarBins);
  float azimuthBinSize;
  float polarCenter, azimuthCenter;
  int numberOfAzimuthBins;

  for(int polarBin=0; polarBin<numberOfPolarBins; ++polarBin)
  {
    polarCenter = polarBinSize * float(polarBin);
    numberOfAzimuthBins = floor(fabs(float(numberOfPolarBins*2) * sin(polarCenter))) + 1;
    azimuthBinSize = (2.f * M_PI) / float(numberOfAzimuthBins);
    for(int azimuthBin=0; azimuthBin<numberOfAzimuthBins; ++azimuthBin)
    {
      NormalType vec;
      azimuthCenter = azimuthBinSize * float(azimuthBin);
      vec[0] = sin(polarCenter) * cos(azimuthCenter);
      vec[1] = sin(polarCenter) * sin(azimuthCenter);
      vec[2] = cos(polarCenter);
      vecs.push_back(vec);
    }
  }
  vecs.push_back(NormalType(0.f, 0.f, -1.f));
  return vecs;
}

